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' Abstract. The mobility problem for suspension of spherical particles immersed in 



an arbitrary flow of a viscous, incompressible fluid is considered in the regime of low 
Reynolds numbers. The scattering series which appears in the mobility problem is 
simplified. The simplification relies on the reduction of the number of types of single- 
particle scattering operators appearing in the scattering series. In our formulation 
\ there is only one type of single-particle scattering operator. 
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1. Introduction 

Various types of suspensions of spherical particles can be found in nature and industry. 
The diversity follows from the facts that systems may be composed of many sorts of 
particles, such as hard-spheres or spherical polymers [1], and many types of forces 
may act between particles [2], [3]. This variety in a structure implies a wide range 
of phenomena exhibited by suspensions. To understand them it is often crucial to 
consider hydrodynamic interactions between particles, that is their mutual influence 
through movement of the surrounding liquid (which is different than the influence 
through direct forces, such as magnetic or van der Waals force). In many physical 
situations examination of the hydrodynamic interactions amounts to the friction 
problem or the mobility problem [1]. In both cases the particles are immersed in a 
flow of the surrounding liquid (ambient flow). In the friction problem the velocities 
of particles are assumed and the hydrodynamic forces on the fluid produced by the 
particles are calculated. This problem appears e.g. in the determination of Stokes 
coefficient for polymer modeled gglomerate of spherical particles [5J El [7] . In the 

mobility problem one determines the velocities of freely moving particles and also the 
hydrodynamic forces acting on the fluid. The particles are assumed to be immersed in 
the ambient flow and the external forces may act on them. Among the situations in 
which the mobility problem appears, the determination of sedimentation coefficient of 
the suspension can be mentioned [El IS]- 

To solve both the friction and the mobility problem one starts with the equations 
which govern the dynamics of the suspension. Here we assume linear, stationary Stokes 
equations for an incompressible fluid with the stick boundary conditions on the surface 
of particles. One of the possible approaches to the Stokes equations is the method 
of successive approximations [10], also known as the reflection method, developed by 
Smoluchowski pj] . It is based on linearity of the Stokes equations. The method consists 
in successive superpositions of the single-particle solutions of the Stokes equations 
chosen to fulfill the boundary conditions with the increasing accuracy. The above 
procedure leads to the solution of the friction problem in a form of the superpositions of 
multiple reflected flows. The resulting structure of the solution is called the scattering 
series. In the series the single-particle operator plays an important role. In case of 
the friction problem the scattering series has a simple form because there is only one 
type of the single-particle operator. The situation is more complicated in the mobility 
problem. Here four types of single-particle operators are considered [12]. Due to the 
difference in the number of single-particle operators statistical physics considerations 
are easier in the case of the friction problem than in the case of the mobility problem. 
For the latter, many formulas of the same structure can be found in the literature 
[13j[Tl[T3[Tg[I2[I3[l8l[T9]. 

The aim of the present article is the reformulation of the mobility problem. The 
reformulation results in the scattering series in which only one type of single-particle 
operator appears. It enables simplification of the statistical physics considerations 
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relevant to the mobility problem. 
2. Governing equations 

The system under consideration consists of N identical hard spherical particles of radius 
a immersed in an incompressible, infinite Newtonian liquid with kinematic viscosity 
77. The inertia of particles and the inertia terms in the incompressible Navier-Stokes 
equations are assumed to be negligible. As a consequence, the fluid is governed by the 
steady Stokes equations 0j. We supplement the Stokes equations with the stick (no-slip) 
boundary conditions on the surface of particles (20]. Next, following idea of Mazur and 
Bedeaux [21 J we extend the Stokes equations inside the particles in the following way: 

Vp (r) — r]Av (r) — f (r) , (la) 
V-«(r) = 0, (lb) 

introducing induced force densities / (r) [22], [23J. Here p (r) and v (r) are the pressure 
field and the velocity field of the suspension. The induced force densities are determined 
pM] by the condition that the flow of the suspension inside the particles reproduces their 
hard-sphere velocity field, 

v(r) = U i (r) = V i + n i x(r-R i ) for \r - R t \ < a, (2) 

where Vi and fij is translational and rotational velocity of the z-th particle which is 
located at the position R4. 

For the case under consideration the force densities / (r) are localized only on the 
surface of particles [IS], [25], that is 

N 

/(r) = $>(r), (3) 

i=i 

with fi (r) localized only on the surface of the z-th particle, 

fi (r) = -cr-n i 5{\r-R i \-a), (4) 

where a is the stress tensor for the fluid [1], rti - a vector normal to the surface of the 
particle i at point r, and 5 (x) - Dirac delta function. 

The above extension of the Stokes equations allows to use the Green function 
method. The method leads to the expression for the flow of suspension in the whole 
space [26] . 

„(r) = J d 3 r'G (r-r')-f(r'), (5) 
where Oseen tensor Gq has the following form [27] 
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3. Friction problem 

In the friction problem the particles are immersed in a liquid in which initially ambient 
flow v (r) is present and their translational VI and rotational velocities fi, are assumed 
to be known. The aim of the friction problem is a determination of the force densities 
fi (r) induced on the surface of particles. 

For a single particle, the friction problem for a particular type of ambient flow was 
solved more than a century ago [10]. It was later generalized for an arbitrary ambient 
flow [21], [21] ■ The solution of the Stokes equations ([1]) in this case has a form of the 
following linear relation [28] 

/! (r) = J d 3 r'Z (R u r, r') ■ {U x (r') - v (r')) . (7) 

Single particle resistance operator Z is localized on the surface of a particle. It means 
that the force density f\ (r) is localized on the particle surface - consistently with the 
equation (@| - and its value depends only on the velocity field U\ (r) — v ( r ) at points 
\r — R\\ = a. The details of the resistance operator Z can be found in the Appendix. 
The equation (j^J) will have the following form in shorthand notation 

/i = Z (l).(L7i-«o), (8) 

in which the integral variables are omitted and the position of a particle is denoted by 
its index. The notation will be extensively used further on. 

It is worth mentioning that the single particle friction problem for a spherical 
shape has been solved not only for a hard sphere with the stick boundary conditions. 
Many other physical situations have also been considered in the literature, e.g. different 
boundary conditions [2H], permeable particles [30], spherical polymers [3Tj . immiscible 
droplets [32], [33] or more complex cases [34], [35]. In general, the relation (jH]) is still 
valid but with modified Zq operator. 

The concept of the induced force densities and linearity of the Stokes equations 
allows to use the single particle friction problem to find the solution of the friction 
problem for suspension. In fact, the i-th particle in suspension is surrounded by a 
flow which is a superposition of the ambient flow vq (r) and the flow induced by other 
particles, Ylij^i I ( r ; r ') fj ( r 0- Applying this modified ambient flow to equation 

(jH]) leads to the following expression 

/i = Z (i). (Ui-v - Y, G °fi) ( 9 ) 

written in the shorthand notation. This is the formula where one can directly implement 
the reflection method by successive iterations. It yields 

TV 

/ i = 5> ij (l...iV)(E/ i -tto), (10) 
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where the friction operator Z^ has the form of the scattering series 

Z iS (1 . . . N) = S iS Z Q (i) - (1 - <y Z (z) G (zj) Z (j) + 

+ £z (i) Go (ifc) Z (fc) Go (fcj) Z (j) + . . . (11) 

k 

The different terms in the equation ( TTTT) correspond to scattering sequences. The prime 
symbol indicates summation over k different than neighboring particle indexes in the 
scattering sequence. 



4. Mobility problem 

In the mobility problem one considers freely moving particles immersed in an ambient 
flow of the fluid v Q (r) and subjected to the action of the external forces. The aim of the 
problem is to calculate the velocity fields of the particles Ui (r) and also the induced 
force densities fi (r) [I] . 

In order to obtain the solution for a suspension we first analyze a case of a single 
particle. Before going into the details, it should be noticed that linearity of the Stokes 
equations implies linear relation between the response of the particle and the source of 
a disturbance. Therefore the velocity field of the particle U (r) and the induced forces 
fi (r) on its surface are linear to the ambient flow Vq (r) and to the external forces 
fext (r), 

C7i (r) = J dr'Mo (R u r, r') f ext (r') + / dr'M < (R u r, r') v (r 1 ) , (12) 

/! (r) = j dr'M > (R u r, r') f ext (r f ) + J dr'M (R u r, r') v (r f ) . (13) 

Since the above statement is based only on the linearity of the governing equations 
it is correct for different particles and boundary conditions mentioned in the previous 
section. The single particle operators Mq, M < , M > and M need to be determined 
for every particular case. In what follows we discuss the hard spheres with the stick 
boundary conditions [TS] . [12] . From the papers cited here, one can easily infer the 
form of M = /x operator which is explicitly given in the Appendix. The remaining 
operators are expressed with the following equations 

M< (1) = Mo (1) Z (1) , (14) 

M>(1) = Z„(1)a*o(1), ( 15 ) 
M (1) = -Z (1) + Z (1) /x (1) Z (1) (16) 

which were written in the shorthand notation. 

Due to linearity of the Stokes equations the above single particle solution of the 
mobility problem can be used to solve the case of suspension. In the suspension, the 
z-th particles is immersed in the flow given by a superposition of the ambient flow v (r) 
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and the flow generated by other particles, Y^jjt% I G o ( r , r ') fj ( r ')- Assuming this 
modified ambient flow in equations (j!2p and ([TBI leads to the expressions for the velocity 
of particles 



U { = M (i) f ext + M < (i) lv + ^2G fA 



;i7) 



and the induced force densities in suspension 



/« = M> (z) / ext + M (i) ( « + ^ Go/, . (18) 
In what follows we rewrite the equations ffTTj) and fllSH in the following concise form 



Mil) 



(19) 



where the response of the z-th particle Sj, single-particle mobility operator M(i), and 
external field i/>o are defined respectively below: 



Ui 
fi 

V 



M(i) - 

whereas G is generalized Oseen tensor G : 



M (i) M<(z) 



M> (z) M (z) 



(20) 

(21) 
(22) 



G 




Go 



(23) 



The method of iterations applied to the equations ([19]) leads to the following 
solution of the mobility problem 



N 



Si (jRi, . . . , Rn) — ^^Tij (Ri, . . . , Rn) V'o? 
where XL- is given by the scattering series as follows 



(24) 



Til (Ri, ■ ■ ■ , R 



N) 



SijM (Ri) + (1 - Sij) M (Ri) GM (Rj) + 

+ M w GM ( R k) gm ( R j) + ■■■ 



(25) 
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It is worth comparing the scattering series given by expression (1251) to the scattering 
series 

<%M (Ri) + (1 - M < (RJ GM {Rj) + J^ M < W G ^ W GM > ^ + ' ' • 

k 

(26) 

which is considered in the literature [TJ]. Notice that here four types of single-particle 
operators M , M<, M>, and M appear. With respect to the number of types of 
single-particle operators, the formulation of the scattering series ( |25l) introduced in the 
present paper is simpler than the series given by expression ( 1261) . 

It is worth mentioning that no approximation was made in the above analysis. 
In particular, the equation ( 1191) is a proper starting point to analyse hydrodynamic 
interactions of particles in close contact. 

5. Discussion 

In the present paper the scattering series for the mobility problem has been reformulated 
which results in the simple form given by the equation ( 1251) . The simplification relies 
on the fact that in expression ( )25|) there is only one type of single-particle operator, M. 
In the formulation hitherto used in the literature [12], in the scattering series there are 
four type of single-particle operators M Q , M<, M>, and M which is showed by the 
expression ( |26l) . 

At first sight the difference between both formulations may not seem to be 
significant. However, the mobility problem plays a crucial role e.g. for calculations 
of transport coefficients of suspensions. In this context statistical physics considerations 
contain many different formulas of the same structure [16], [12], [13], [TJ]. A forcible 
example of cumbersomeness following from a lack of a simple formulation of the mobility 
problem is the Beenakker-Mazur method [17] , [18] , [19] which is used to calculate short- 
time dynamic properties of suspensions. A multitude of expressions occurring in cited 
articles obscures the essence of the method which, on the other hand, is the most 
comprehensive analytical scheme available so far [36]. With the aim of the reformulated 
scattering sequence the Beenakker-Mazur method can be presented in a simple form 
which we are going to show in another paper. The simple formulation of the mobility 
problem also allows to carry out advanced analysis of scattering series in a clear and 
simple way. In this context the formulation will also be used in our subsequent work on 
macroscopic characteristics of suspensions. 
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Appendix A. Single particle operators and multipole picture 

Here we give the explicit form of Zq and operators. Basing on the reference [25] the 
Z operator can be expressed by the following formula 

oo {' I 2 

Z (R,r,r')= E E E 5 *( r - R ) w L*( r - R ) 

l,l'=l m'=-V m=-l a,a'=0 

x [Zo (•R)] fcnff> ,/ m v S a (r' - R) wr m , a , (r' - R) , (A.l) 

where [Zo (R)] lm(Tl , ml(T , stands for the multipole matrix with indexes I = 1, . . . , oo; 
m = —I, cr = 0, 1, 2. Its matrix elements are explicitly given e.g. in the reference 
[37] . A set of multipole functions v^ ma (r) and w^ ma (r) is defined e.g. in the references 
[38] or [25]. Every solution of the homogeneous Stokes equations may be expressed as a 
combinations of the multipole functions vt. (r) 

oo I 2 

"oM = EEE N ( R )\ lm , (r-R), (A.2) 

1=1 m=—l a=0 

whereas the induced surface force fi (r) as combination of multipole functions w^ ma (r): 

oo I 2 

/«W = EEE Mima S * ( r ~ R ) ™Lr ~ R ) ■ (A.3) 

1=1 m=—l a=0 

The multipole functions w^ ma (r) are defined as orthonormal to (r) functions. This 
orthonormality is expressed in the following way 

( S awf ma \v$ m , al ) = 5 w 5 mm ,5 call (A.4) 

with the Dirac notation [39] for scalar product of two vector fields A (r) and B (r): 

(A\B) = J d 3 rA* (r)-B{r), (A.5) 

and the scalar function 5 a (r) of the form 

5 a (r) = a' 1 6 (|r| - a) 

which confines integration area to the sphere of the radius a. 
Operator /n is given by the expression 

oo V I 2 

MH.r.r'H E E E E Q *( r - R ) v L,(r- R ) 

l,l'=l m'=-V m=-l a,a'=0 

x [Mo («) - - «) v /w - «) ( A -6) 

with multipole matrix [P'o(R)]i ma i< m t a i explicitly given in the reference [37]. Here 
G a (r — -R) is a characteristic function of the particle at position R: it equals whenever 
r points outside the particle, and equals 1 otherwise. 
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Finally, we represent the equations (j!7p and (1181) in the multipole expansion 
formalism. To pass on to the multipole picture we put the expressions flA.lj) and flA.6j) 
into these equations and multiply them by {w^ mcr (i) 8 a (i) | from the left side. A simple 
algebra yields 



Ui = no (i) fext (i) + Ho («) Z (i) v (i) + ^ G ° (#) fi 



fi = Z (i) no (i) f ext (i) - Zo (i) ^ (i) + ^ G ° (#) f^j 
where the multipole vector of the ambient flow at point R 

[M4» = K-( fl )Mfl)l%), 

multipole velocity field Ui for the z-th particle 
induced surface force multipole fi for the i-th particle 

[Mima = ( V La (*) > 

and external force multipole field / ext ( J2) at point i? 



(A.7a) 
(A. 7b) 

(A.8) 
(A.9) 
(A.10) 
(A.ll) 



In the above formulas | A (i)) or \A (R)) denote vector fields A(r — R4) or A(r — R) 
respectively. Moreover matrix Go {R, R') is defined with the formula 

[Go (R, R% m a,l'm>*> = {™tna («) *a («) | ^0 |tO? mV (#) 5 a (#)> . (A.12) 

Its matrix elements may be found in the references [37] , jlQ] . It is worth mentioning that 
[G (R, R')]i mcT vm'a 1 depends on the difference of positions R—R' and for nonoverlapping 
configurations, i.e. \R— R'\ > 2a, it scales as 1/\R—R 



nW '+&+(?' '— 1 



equations ( 1A.7j) we also used the following definition [16] 

Zo (i) = Z (i) - (i) /-to (») (i) 



To obtain 



(A.13) 



In the multipole formalism the integral equations (jA77j) may be easily reformulated 
into equation 



Si = M (Ri) I tfj (Ri) + J2 G R ^ s i J 



(A.14) 



in the same way, as equations ( IT71) and ( ITS!) were transformed into the equation ( TK) 
Moreover definitions of Sj, an d are similar to definitions fl20j) . ( 12 ip and f l22|) : 



(A.15) 
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4>o(R) 



fext (R) 

v (R) 



M(i) 



fio (i) /J (i) Z (i) 
Z (i) A*o (i) -^o (i) 



and the Green function G in extended multipole space has the following form: 



G 





Go 



(A.16) 
(A.17) 

(A.18) 
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